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Abstract 

Short-  and  long-time  solutions  for  material  balance  equations  for  porous  electrodes  in  both  the  solid  and  the  solution  phase  is  presented. 
These  solutions  represent  the  concentration  profile  of  lithium-ions  in  the  solid  and  solution  phases  of  the  positive  electrode  of  a  lithium-ion 
cell,  and  are  obtained  using  the  method  of  Laplace  transform  for  short-  and  long-time  intervals. 
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1.  Introduction 

This  paper  presents  an  analytical  solution  for  the  material 
balance  equation  that  governs  the  concentration  distribu¬ 
tions  in  the  positive  electrode  (cathode)  of  a  lithium-ion  cell 
operating  under  galvanostatic  conditions.  It  is  an  extension 
to  the  work  reported  by  Atlung  et  al.  [1]  who  have  derived 
analytically  a  short-time  solution  for  the  lithium-ion  con¬ 
centration.  These  workers  [2]  have  also  obtained  a  short- 
time  analytical  expression  for  the  discharge  curve  in  a 
system  using  composite  insertion  electrodes. 

Recently,  Subramaniam  and  White  [3]  presented  an  exten¬ 
sion  to  the  method  of  separation  of  variables  for  solving  the 
equations  that  control  the  concentration  distributions  in  solid 
electrodes.  In  this  work,  we  have  solved  the  material  balance 
equation  by  the  method  of  Laplace  transform.  The  material 
balance  equation  takes  into  account  the  migration  term,  but 
for  a  constant  lithium-ion  transference  number.  The  solutions 
obtained  are  for  short-  and  long-time  intervals. 

2.  Description  of  system 

The  lithium-ion  cell  that  is  being  considered  in  this  work 
has  the  configuration  LiMn204 1 separator  I  graphite  with  a 
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liquid  electrolyte.  A  schematic  representation  of  a  clo¬ 
sely-packed  cathode  employing  LiMn204  active  material 
is  shown  in  Fig.  1.  The  particles  of  active  material  are 
considered  to  be  spherical  and  have  high  porosity.  The 
electrodes  are  composites  that  implies  that  the  active  mate¬ 
rials  are  also  mixed  with  carbon  black  to  increase  the 
electronic  conductivity.  This  model  uses  the  porous  elec¬ 
trode  theory  [4]  and  the  concentrated  solution  theory  [5]  to 
treat  the  composite  electrodes.  During  the  discharge  process, 
lithium-ions  deintercalate  from  the  negative  electrode 
(anode),  diffuse  through  the  separator,  and  intercalate  into 
the  cathode. 

3.  Solution  phase  of  cathode 

In  order  to  compare  our  results  with  others,  we  have 
written  the  equation  in  accordance  with  that  of  Fuller  et  al. 
[6],  where  the  governing  equation  for  the  electrolyte  con¬ 
centration  is  given  as: 

+(x,,)  =  s^£e(-«,.)  +  ++  <D 

where  0(x,  t)  is  the  material  balance  of  salt  in  the  solution 
phase.  Here,  D  is  the  diffusion  coefficient  of  the  lithium-ion 
in  the  electrolyte,  e  the  volume  fraction  of  the  solid,  t+  the 
transference  number,  a  the  specific  interfacial  area,  and  jn  is 
the  pore  wall  flux  of  lithium-ion.  Similar  to  Doyle  and 
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^ -  anode - ►  separator  ^ -  cathode - ► 


x=()  x=8s  x=8c 

Fig.  1.  Schematic  representation  of  lithium-ion  battery  used  in  this  model. 


Newman  [7],  the  following  dimensionless  parameters  are 
introduced: 


©  x  Sc  Dt 


Eq.  (1)  when  written  in  its  dimensionless  parameters 
becomes: 


dc 

dz 


+  J 


subject  to  the  boundary  and  initial  conditions 

dc 

--  =  0  at  y  =  1  +  r 
dy 

dc  Jr 

=  al),=  1 

and 


(2) 

(3) 

(4) 


c(y,0)  =  l  (5) 

The  term  ajn  in  Eq.  (1)  represents  the  rate  of  transfer  of 
lithium-ions  from  the  electrolyte  solution  inside  the  pores  to 
the  solid,  and  is  related  to  the  divergence  of  the  current  flow 
in  the  electrolyte  through  Faraday’s  Law  [8],  where  an 
average  of  jn  [7]  is  taken  to  be: 


-I 


Jn 


aFL 


or  J  = 


+ 


zip  -  4)Z 

FDcosr 


(6) 


On  application  of  the  Laplace  transformation  with  respect  to 
dimensionless  time  z,  we  get: 


dc 

dz 


d2c 

dy2 


■OO 


dz  F  J 


-ST 


dz 


Integrating  both  sides  of  Eq.  (7)  and  interchanging  the  orders 
of  differentiation  and  integration  of  the  first  term  on  the 
right-hand  side,  we  obtain: 


d2c  s  _  1 

dy2  D  81/2 


which  gives  the  solution  of: 


/ 


c(y,  s)  =  A  cosh(gy)  +  B  sinh(gy)  +  1  -j-  - 


(9) 


where  q  =  \J s/e,1/2  and  the  boundary  conditions  are: 
dc 


d y 

dc 


0  at  y  =  1  +  r 
Jr 

at  y  =  1 


(10) 


dy  e}/2s 

Applying  the  boundary  conditions  (10)  and  (11)  gives: 
Jr  coshg(l-j-r) 


(11) 


A  =  - 


e}/2sq  sinh(gr) 


(12) 


and 


B  = 


Jr  sinh^(l-f-r) 


rl/2sq  sinh(gr) 

Eq.  (9)  now  becomes: 

,  .  \  J  Jr  /coshg(l  +  r  —  y) 

ay,*)  =  -  +  -?- 

s  sJ 


(13) 


>2  e1/4  y  s2/2  sinh(gr) 


(14) 


Expanding  the  trigonometric  term  in  Eq.  (14)  in  its  expo¬ 
nential  form,  Eq.  (14)  can  be  rewritten  as 

1  /  Jr  f  Q-q{J-y+2r)  e-?(y-i)' 

s)  =  -  + -5  -  75J 


5,3/2(l  -  e_2^r) 

1  /  Jr  1  q(\-y+2r\n+\\) 

S  S2  81/4  S2/2  ^ 


+  e 


-q(y—  l+2m) 


n= 0 


(15) 


Using  the  Laplace  transformation  Table  [9],  we  get: 


.  f coshg(l  +  r-y) 
\  si!2  sinhfr//') 


OO 


=  E)  (4t) 


1/2 


n= 0 

+ 


i  erfc 


(1  -y  +  2r[n+  1])' 


2  V  s1/2- 


■•erf Tly-l+j"-)) 

V  2V  81/2!  ) 


(16) 


where: 


1 


i  erfc(v)  =  — — -  e  *  —  verfc(v), 
—  V - 


71 


(17) 


Hence,  the  final  solution  of  Eq.  (15)  is  given  by: 

2Jrz1!2  F^F  ( _  -  /I  —  y  +  2r[n  F  1] 


c(y,  t)  =  1  +  Jz  — 


■1/  4 


(  i  erfc 


n= 0 


2e1/4i1/2 


+ 


i  erfc  ( (y  ~  1  +  2m))  ) 
V  2s1/4Tl/2  )) 


(18) 


It  can  be  observed  that  Eq.  (18)  is  valid  only  for  a  short  time. 
In  fact,  for  large  values  of  s ,  Eq.  (15)  will  give  the  con¬ 
centration  profile  as: 


c(y,  r)  =  I+Jt 


(19) 
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and  the  concentration  increases  linearly  with  time  for  a 
particular  value  of  y. 

In  order  to  obtain  the  long-time  solution,  the  trigono¬ 
metric  term  of  Eq.  (14)  is  expanded  as  follows: 


cosh  z=  1  + 


4z: 


7T 


sinhz  =  z  1  + 


7 v 


1  + 


1  + 


4? 

32n2 

227l2 


1  + 


1  + 


4? 

52n2 

327l2 


(20) 


(21) 


Defining: 

f{s)  cosh  q(  1  +  r  —  y) 
g(s )  s 3/2  sinh(gr) 

coshg(l  +  r  —  y) 

s2{  1  +  (q2r2 / 7T2) ) ( 1  +  (q2r2  / 22n2)) . . . 


(22) 


where  f(s)  and  g(s )  are  polynomials  in  s  which  have  no 
common  factor  and  the  degree  of f(s)  is  lower  than  the  degree 
of  gCs). 

Clearly,  s  =  0  is  a  root  for  Eq.  (22)  and  the  other  root  is 
given  by  q  =  ( nni)/r  or  equivalently  s  =  (n2n2sl^2)/r2 
where  n  =  1,  2,  3,  . . . 

By  partial  fraction,  Eq.  (22)  becomes: 


coshg(l  +  r  —  y) 
s2/2  sinh(gr) 

n  =  1 , 2,  3, . . . 
This  implies  that: 

cos  q(  1  +  r  —  y) 


oo 


A  2  ^ — > 

^4-  —  +  Vtt 

s2  s  '  ( 1 


B 


m 


m= 


'  (1  +  (q2r2 / n2n2))  ’ 


(23) 


oo 


Msft  (  1  + 


n=  1 
oo 


2  2' 
ti2tc2 


+  A2J^[ (  1  4~ 


n= 1 


<?2r2 

ti2tc2 


00 


00 


+ s2yy?m  n  ( 1  + 

m=  1  n=l,nfm 


2  2 
n2TC2 


(24) 


Taking  the  limit  as  s  — >  0,  we  get  A2  =  1. 

Eq.  (24)  can  now  be  rewritten  as: 

coshg(l  +  r  —  y) 

=  SMl(|r)+/1|S„2siiih(<;r)+j2g^  pj 


l 


772=1  n=l,nfm 


tends  to  zero,  and  applying  L’Hopital’s  rule,  we  get: 


Ai  = 


1  /  (1  +r-y)‘ 


1  +  r 


4/4 


2  r 

U72 


Letting  q  =  ( kni)/r  and  substituting  it  into  Eq.  (25)  gives: 


cosh 


0  +  r-y)f) 


:l/2k27ix 


00 


Bm  ]^[ 

n=l,nfk 


1  - 1  - 

n 


£+1 


or 


kn  ,  sk4n4  (— 1) 

cos  —  (1  +  r  -  y)  =  — — Bm - - 


Eq.  (27)  implies  that  the  values  for  Bm  are: 
(— l)A+12r4cos(k7i/r)(l  +  r  —  y) 


(27) 


Bm 


&k47i4 


(28) 


Hence,  the  partial  fraction  in  Eq.  (23)  becomes: 
coshg(l  +  r  —  y) 


s3/2  sinh(gr) 


1  1 

—  ~2  +  “ 

sz  s 


2 81/4  ( (1  +  r  —  y )2  r 


00 

+  E 


1  +  r  \  2s1/2  e3/4 

(— l)n+12r4  cos(mc/r)(l  +  r  —  y) 


(29) 


n_x  m4n4(\  +  (q2r2  /  n2n2)) 

Eq.  (29)  is  then  inserted  into  Eq.  (14).  Upon  taking  the 
inverse  of  the  Laplace  transform  of  the  substituted  Eq.  (14), 
we  get: 

c(y,  t) 


=  1  +  /T  1  - 


2 Jr  ( (1  H-  r  —  y)‘ 


Jr 


81/47  1  4-  r 

2r\  ,2 


2e‘/2 


;3/4 


7l2g3/4 


E^  ^  exp(— A2  81/2t)  cos(/l2  (1+r— y)) 

n=\ 


(30) 

where  2m  =  (mn)/r.  This  equation  is  valid  for  a  long-time 
interval  and  is  related  to  the  discharge/charge  coefficient. 


4.  Solid  phase  of  electrode 

In  the  solid  phase,  the  porous  electrode  solid  particles  are 
considered  to  be  spherical  with  a  radius  R.  The  concentration 
of  lithium-ions  inside  a  solid  sphere  of  active  material,  cs ,  for 
a  definite  radius  r  and  time  t  is  given  by: 

rd2c 

'  "  '  (3D 


dc 


dt 


=  D, 


_ >5  2  dc, 

<9r2  r  <9r 


'2r2\ 

As  in  [10],  Eq.  (31)  is 

subjected  to  the  following  initial  and 

2n2) 

boundary  conditions: 

CS(V,0)  =  Co 

(32) 

(25) 

l  as  .v 

dcs(0,t)  _  0 

dr 

(33) 

(26) 

n  dcs(R,  t)  _  . 

0  Jn 

or 

where  c0  is  the  initial 

(34) 

concentration  of  lithium-ion  in  the 

solid  particle,  jn  the  average  pore  wall  flux  of  lithium  at  the 
surface  of  the  particle  [11],  and  Ds  the  diffusion  coefficient 
of  the  lithium-ion  inside  the  particle.  The  negative  sign 
indicates  that  the  pore  wall  flux  across  the  solidlliquid 
interface  is  in  the  opposite  direction  to  the  lithium-ion 
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diffusion  into  the  surface  of  the  insertion  material.  The 
radius  of  the  particle  is  traversed  by  the  lithium-ion  when 
r  =  R.  Atlung  et  al.  [1]  and  Subramaniam  and  White  [3] 
have  set  a  zero  initial  condition  in  their  work. 

Let  y  =  r-c,  then  the  Laplace  transform  of  the  linearized 
Eq.  (31)  is  given  by: 


dr2 


s  _  c0 

—  y  = - r 

D/  Ds 


(35) 


Eq.  (35)  was  integrated  with  respect  to  r  twice,  which  is  why 
Eq.  (35)  becomes 


From  Eq.  (42),  it  can  be  observed  that  the  Laplace  transform 
can  only  be  inverted  for  small  time  intervals.  Upon  inversion 
[9]  for  n  =  1  and  k  =  1,  we  get: 


Cs(P,0 


Co  + 


Rjn 


Dsp 


—  exp(—  1  +  p  +  t)  erfc  (  — 


i  +  /A 

2^) 


-ferfc 


Af) +exp(- 


—  1  —  p+T)  erfc  (  —  vt  + 


1  + 


£\ 


-  erfc  ( l-Af) 


vvU. 


2\A  ) 

(43) 


y{r,s)  =  B{tqr  -  t~qr)  +  —  r 

s 

or 


(36) 


where  the  two  dimensionless  variables  are  defined  as: 


and 


T  = 


Dst 
R 2 


rc(r,s)  =  2Bsinh(qr )  +  —  r  (37) 

s 

where  q  =  s/Ds. 

Here,  B  is  one  of  the  two  constants  of  integration  obtained 
upon  performing  the  integration  twice  on  Eq.  (35).  It  can  be 
easily  shown  that  the  other  constant  of  integration  is  —B. 

On  differentiating  Eq.  (36)  and  applying  the  boundary 
condition  (34),  we  have: 


2  B 


sDs&r[(l  —  qr)  —  (1  +  qR)  e-2^] 


(38) 


Hence,  Eq.  (37)  now  becomes: 


Eq.  (43)  gives  a  relationship  between  the  lithium-ion  con¬ 
centrations  in  the  solid  phase  with  discharge  current  only  for 
short-time  intervals. 

On  taking  the  Laplace  transform  to  Eq.  (31),  we  have 
arrived  at  Eq.  (37).  Upon  applying  the  boundary  condition 
(34),  we  obtained  Eq.  (38),  which  we  now  rewrite  as 


2  B  = 


JnR2 _ 1 _ 

sDs  sinh (qR)  —  qRcosh(qR) 


(44) 


Hence,  Eq.  (37)  becomes: 


'  =  Co  jnR2 _ sinh(gr) _ 

s  rDs  s[sinh(qR)  —  qRcosh(qR)} 


(45) 


—  /  \  Co 

c(r,s) - b 
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s  rsDs{  1  —  qR) 


(e 


~q(R~r)  _  ~-q(R+r) 
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c{r,s) 
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1  _  1  +_gg  e-2 qR 


“I  -1 
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1  -  qR 
JnR2 


°°  /\+qR\  __2  nqR 


c(r,s) 


s  rsDs(  1  —  qR) 

Co  ,  jnR 2 


(e 


-q(R-r)  _  „-q{R+r 


s  +  rDs 


(e- q{R~r )  _  Q~q{R+r)^ 


oo  n 


X 


n 


^ _ y)k2n~k  Q—^nqR 


n= 0  k=0 


k  )  s(  1  —  qR) 


n+l—k 


,  ,  Co,  jnR2  [  -  e~*R+r 

c(r,s)  = - h 


s  rD , 


+ 


s(  1  —  qR) 

1 


(39) 


(40) 


(41) 


s(l-qR)  s(l-qR) 


(e 


-q(3R-r)  _  -q(3R+r) 


) 


+ 


4  1 

+ 


s(l  —  qR)3  s(l  —  qR)2  s(l  —  qR) 


x 


^Q-q{5R-r)  _  Q-q(5R+r ))  _ 


(42) 


We  must  now  consider  the  inverse  transform  of  the  second 
term  on  the  right-hand  side  of  Eq.  (45).  Making  use  of 
Maclaurin  series,  we  have: 


-2  R3  y^(n+l)RZns 


2 n 


sinh  (qR)  —  qRcosh(qR)  — 


A 


n= 0 


(2n  +  3)!  D 


(46) 


Hence,  the  partial  fraction  of  this  second  term  is  given  by: 
sinh(gr) 

sinh  (qR)  —  qRcosh(qR)] 

qr[  1  +  (M2/ 3!)  +  ((<?r)4/ 5!)  +  •  •  •] 
{—2R}/Ds)qs2  £((«  +  \)R2n)/{{2n  +  3)!Ds"s«) 


sinh(gr) 


sinh(qR)  —  qRcosh(qR)}  s 
where: 


G(s) 


(48) 


sinh(gr) 


Q 


Aq  G(s)  +  A\s  G(s )  +  s2  F(s) 


(49) 


and 


—2 R3  (n-\-l)R2n  n  sinh(qR)  —  qRcosh(qR) 

(>=-5rh<^ w' = - * - 


(50) 
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If  Xn  are  the  roots  of  G(s)  =  0,  then  Xn  —  —Ds(x2/R2,  where 
ccn  are  the  non-zero  roots  of  tan  6  =  0.  From  Eqs.  (49)  and 
(50),  we  have: 


(51) 


A 


l 


3 r  ( R 2  r2\ 

¥  \To  -  ~6j 


(52) 


and 


no 

G(s) 


OO 


*E 

n=  1 


sin((«„/r)/j?)  1 

a2  sin  an  s  —  A„ 


(53) 


where  values  are  simple  zeros  of  an  nth  order  polynomial 
G(s). 

Hence,  the  Laplace  transform  of  Eq.  (45)  can  be  written 
as: 


cs(r,t)  =  co  ~ 


J„n  3  Dst  1  /r\  2  3  JR 
Dr  1  R2  +  2  \RJ  10  V  r 


x 


OO 


E 


sin(mn/R) 


n=1  a-sina„ 


exp(-alDst/R2) 


(54) 


where  the  ocn  values  are  the  positive  roots  of  ocn  cot  ocn  =  1 . 
Eq.  (54)  is  valid  for  all  t.  As  t  tends  to  infinity,  the  summation 
tends  to  zero  and  Eq.  (54)  is  simplified  to: 


cs(r,t)  =  c0  - 


j„R2  (3Dst  1  (ry  3 
1  R2  '  2  \r)  10 


(55) 


5.  Results  and  discussion 

The  profiles  of  the  concentration  of  lithium-ions  during  a 
galvanostatic  discharge  in  the  solution  phase  of  the  cathode 


Fig.  3.  Concentration  profile,  for  long-time  interval  in  solution  phase 
across  cathode  during  galvanostatic  discharge  of  1.39  mAcm-2. 


are  plotted  in  Figs.  2-4.  Fig.  2  is  plotted  using  Eq.  (18),  and 
Figs.  3  and  4  using  Eq.  (30).  The  solid  volume  fraction  e  is 
taken  to  be  0.524  [7]  with  a  total  thickness  of  cathode  of 
v  =  0.174  x  10  4  cm  [7],  and  the  initial  concentration  is 

_ o 

Co  =  3.9  mol  cm  .  From  Fig.  2,  the  concentration  of 
lithium-ions  in  the  solution  phase  is  lowest  at  the  back  of 
the  cathode.  This  is  reasonable  because  during  full  charge,  it 
is  expected  that  the  lithium-ions  concentration  in  the  cathode 
will  be  almost  zero,  and  during  discharge  the  lithium-ions 
fill  up  the  back  portion  of  the  cathode  last.  The  concentration 
of  lithium-ions  at  every  point  in  the  cathode  increases  with 
the  increase  in  discharge  current  for  a  particular  time. 


Fig.  2.  Concentration  profile  for  short-time  interval  in  solution  phase 
across  the  cathode  during  galvanostatic  discharge  of  I  =  1.0  and 
1.39  mA  cm-2. 


Fig.  4.  Concentration  profile  for  solution  phase  across  cathode  during 
discharge  of  /  =  0.5,  1.0  and  1.39  mA  cm-2. 
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Fig.  5.  Concentration  profile  in  solid  phase  across  cathode  during 
galvanostatic  discharge  for  short-time  interval. 


The  filing  up  of  the  cathode  for  a  constant  discharge 
current  at  various  times  is  shown  in  Fig.  3.  For  t  =  3000  s, 
the  cell  is  not  fully  charge  since  the  discharge  current 

_ o 

(/  =  1.39  mA  cm  )  sends  the  lithium-ions  until  y  =  3.12 
into  the  bulk  of  the  cathode.  As  time  increases,  more  and 
more  lithium-ions  enter  the  cathode  and  reach  the  back  of 
the  cathode.  As  time  increases  further  to  5400  and  7200  s, 
the  concentration  of  lithium-ions  increases. 

Plots  of  the  discharge  curves  in  the  solution  phase  for 
various  discharge  currents  at  fixed  time  are  given  in  Fig.  4. 
Again,  the  concentration  is  minimum  at  the  back  of  the 
cathode,  that  is  at  y  =  1  +  r.  The  profile  agrees  with  the 
results  obtained  by  Doyle  and  Newman  [7],  but  the  inter¬ 
section  point  of  the  graph  occurs  near  y  =  3  compared  with 
that  of  Doyle  and  Newman,  where  the  intersection  point 
occurs  at  y  =  2.  This  is  the  effect  of  taking  the  boundary 
condition  at  the  separator  I  cathode  interface  as  Jr/e  ,  where 
J  is  the  dimensionless  migration  term.  Doyle  and  Newman 
have  assumed  a  steady-state  situation  when  plotting  this 
graph. 

Figs.  5  and  6  are  obtained  by  plotting  Eqs.  (43)  and  (54). 
Fig.  5  gives  the  plot  of  concentration  versus  distance  for  a 
short-time  interval  in  the  solid  phase  of  the  cathode  at  a 

_ o 

constant  discharge  current  of  1.39  mA  cm  .  This  plot 
shows  a  hyperbolical  increase  in  concentration  as  the 
lithium-ions  move  from  the  center  to  the  surface  of  the 
electrode  particle.  The  concentration  profiles  inside  solid 


Fig.  6.  Concentration  profile  in  solid  phase  across  cathode  during 
galvanostatic  discharge  for  long-time  interval. 


particles  for  a  large  time  interval  at  a  constant  discharge 
current  of  1 .39  mA  cm  ~  are  given  in  Fig.  6.  Similar  to  Fig.  5, 
the  plot  suggests  that,  during  discharge,  lithium-ions  inside  a 
solid  particle  move  away  from  the  center  of  the  spherical 
particle.  This  is  due  to  the  fact  that  lithium-ions  are  moving 
towards  the  surface  of  the  particle  in  order  to  be  dispersed  into 
the  electrolyte  in  the  cathode  region.  The  concentration  of 
lithium-ions  increases  in  a  consistent  rate  and  becomes  more 
rapid  as  time  increases.  For  a  very  small  time  interval,  the  rise 
in  lithium-ion  concentration  is  rather  slow,  as  shown  in  Fig.  5 
as  compared  with  the  increase  during  a  longer  time  interval. 
This  is  due  to  the  fact  that  a  very  small  amount  of  lithium-ions 
have  been  dispersed  in  this  short  period. 

Even  though  there  has  been  a  lot  of  work  on  mathematical 
modeling  of  lithium-ion  batteries,  most  of  the  work  has  been 
done  numerically  [5,6,8,11-15].  Very  few  analytical  studies 
are  available  in  the  literature.  Only  recently  was  the  material 
balance  equation  for  the  solution  phase  solved  by  involving 
the  migration  term  analytically,  but  this  was  achieved  by  the 
method  of  separation  of  variables  [3,7]  and  not  by  Laplace 
transform  as  in  the  present  work.  Also,  the  solid  phase  has 
been  solved  by  Laplace  transform  for  the  initial  condition 
where  the  lithium-ion  concentration  in  the  sphere  is  zero. 
The  present  work  considers  a  constant  lithium-ion  concen¬ 
tration  in  the  solid  particle. 

The  following  table  differentiates  the  work  done  by  other 
authors  with  that  performed  in  the  present  work. 


Atlung,  West  and  Jacobsen  [1] 


Solved  analytically  the  following  equation  for  the  solid  phase: 

dc 


dt 


d2c  2  dc 
dx2  +  xdx 


0  <  x  <  r 


The  initial  and  boundary  conditions  were: 
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Doyle  and  Newman  [7] 


Doyle  and  Newman  [7] 


Subramaniam  and  White  [3] 


Subramaniam  and  White  [3] 


Present  authors 


(i)  for  t  =  0,  c  =  0, 

(ii)  *  - 


0,  47 ix2  —  =  0, 
ox 


.....  dc  j0 

(ill)  x  =  r,  —  =  —  — . 

ox  D 

Method  of  solving:  Laplace  transform  solution  obtained  for  short-time  only. 
Solved  analytically  the  following  equation  for  solution  phase: 


de 


2  =  e>/2 


o2e 


H-  J 


Ot  "  dy 2 

The  initial  and  boundary  conditions  were: 
(i)  62  =  0  at  t  =  0, 

002 


(ii) 


dy 


0  at  y  =  1  H-  r, 


(iii)  £3/2 


00  2  00 1 


at  y  =  1, 


<9y  <9y 

(iv)  02  =  0!  at  y  =  1, 

Method  of  solving:  separation  of  variables  Solution  obtained  for  long-time  only. 

Solved  the  following  equation  for  solid  phase: 

0cs  \02cs  2  0cs 

_ _  =  n  _ _  j _ _ 

Ot  s  Or 2  r  Or 

Subjected  to  initial  and  boundary  conditions: 

(i)  cs  =  c°s  at  t  =  0, 

Or 

(ii)  =  0  at  r  =  0, 


Or 

(iii)  -A 


0cs 


at  r  = 


<2^+ 

Method  of  solving:  separation  of  variables  Solution  obtained  for  short-  and  long¬ 
time  intervals. 

Solved  analytically  the  equation  in  the  solution  phase  given  by: 


00 


2  =  e1/2 


dl9: 


+  J 


0t  0y2 

The  initial  and  boundary  conditions  were: 
(i)  02  =  1  at  t  =  0, 

00  2 


(ii) 


Oy 


0  at  y  =  1  +  r, 


(iii)  e3/2 


00  2  00 1 


at  y  =  1, 


<9y  Oy 
(iv)  02  =  0\  at  y  =  1 . 

Method  of  solving:  separation  of  variables  solution  obtained  for  short-  and  long¬ 
time  intervals. 

Solved  diffusion  in  spherical  electrode  particle  given  by  equation: 

Oc  ( 02c  2  Oc 

Ot  ydr2  ^  r  Or 

Subjected  to  initial  and  boundary  conditions: 

(i)  c(r ,0)  =  0, 

(ii)  |  (0,  t)  =  0, 

(iii  )%At)  =  Z. 

Method  of  solving:  separation  of  variables  solution  obtained  for  short-  and  long¬ 
time  intervals. 

Solved  analytically  concentration  of  lithium-ions  in  solution  phase  of  a  cathode 
given  by  equation: 


Oc 


02c 


442 


S.A.  Hashim  Ali  et  al./  Journal  of  Power  Sources  112  (2002)  435-442 


Present  authors 


Subjected  to  the  boundary  and  initial  conditions: 

dc 

(l)  at  —  =  0  at  y  =  1  +  r, 
ay 

dc  Jr 

<">5  =  ?7 5a'J'=1- 

(iii)  c(y,  0)  =  1 

Method  of  solving:  Laplace  transform  Solution  obtained  for  short-  and  long-time 
intervals. 

Solved  equation  for  concentration  of  lithium-ions  in  solid  phase  of  cathode  is  given 
by 

dcs  ( d2cs  2  dcs 

Hl  =  D\-dA  +  -rlfr 

Subjected  to  the  initial  and  boundary  conditions: 

(i)  cs(r,  0)  =  c0, 

(ii) ^  (o,o  =  o. 


dr 


dc 


(iii)  -Ds^(R,t)  =j 

Method  of  solving:  Laplace  transform  solution  obtained  for  short-  and  long-time 
intervals. 


6.  Conclusions 

This  paper  presents  the  analytical  solution  for  the  concen¬ 
tration  profile  of  lithium-ions  in  the  solid  and  the  solution 
phases  of  the  positive  electrode  of  a  lithium-ion  cell.  The 
method  of  Laplace  transform  has  been  used  extensively  in 
solving  the  material  balance  equation  in  both  phases  of  the 
cathode  for  both  short-  and  long-time  intervals,  and  thus, 
extends  the  work  of  Atlung  et  al.  [1].  In  solving  for  the 
concentration  profile  in  the  spherical  particle,  the  initial  con¬ 
dition  is  taken  to  be  c0  instead  of  a  zero  initial  condition  as  used 
by  both  Carslaw  and  Jaeger  [16]  and  Atlung  et  al.  [1].  Again, 
we  have  extended  the  work  of  Atlung  et  al.  [1,2]  by  deriving  the 
analytical  solution  for  both  short-  and  long-time  intervals.  The 
explicit  solutions  for  short-  and  long-time  intervals  of  both 
phases  are  presented  together  with  their  graphs. 
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